NSF-KITP-08-94 
MADPH-08-1512 

June 2008 



Spin Dependence of Dark Matter Scattering 

Vernon Barger a ' b J*] Wai-Yee Keung c J3 and Gabe Shaughnessy^ 
"■Department of Physics, University of Wisconsin, 
1150 University Avenue, Madison, Wisconsin 53706 USA 

h Kavli Institute for Theoretical Physics, 
University of California, Santa Barbara, CA 93106 USA 
c Physics Department, University of Illinois at Chicago, Illinois 60607-7059 USA 



Abstract 

New experiments designed to discover a weakly interacting dark matter (DM) particle via spin 
dependent scattering can distinguish models of electroweak symmetry breaking. The plane of spin 
dependent versus spin independent DM scattering cross sections is a powerful model diagnostic. 
We detail representative predictions of mSUGRA, singlet extended SM and MSSM, a new Dirac 
neutrino, Littlest Higgs with T-parity (LHT) and Minimal Universal Extra Dimensions (rnUED) 
models. Of these models, the nMSSM has the largest spin dependent (SD) cross section. It has 
a very light neutralino which would give lower energy nuclear recoils. The Focus Point region of 
mSUGRA, rnUED and the right handed neutrino also predict a very large SD cross section and 
predict a large signal of high energy neutrinos in the IceCube experiment from annihilations of 
dark matter in the Sun. We also describe a model independent treatment of the scattering of DM 
particles of different intrinsic spins. 
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I. INTRODUCTION 



The nature of the dark matter (DM) that comprises about 23% of the energy density 
of the Universe is unknown and a large experimental effort is devoted to its explication. 
The premise is that dark matter is a stable elementary particle. The large scale structure 
of the Universe excludes hot dark matter as the primary DM component. There are a 
number cold and warm dark matter particle candidates, many of which are well motivated 

by theoretical models of TeV scale physics. For recent reviews see Refs. 0, y, y, y, 3. 

Some DM candidates (e.g. axino, heavy gravitino, keV sterile neutrino) may have only 
cosmological consequences. The axion, which can solve the strong CP problem, may be 
detectable from their conversion by photons in a strong magnetic field. The DM particle of 
models with a conserved discrete symmetry (e.g. the stable particle of supersymmetry with 
a conserved R-parity) can have rich consequences for collider experiments, direct detection 
experiments (the elastic scattering of DM particles on nuclei), and astrophysical experiments 
(via the detection of gamma rays, positrons, antiprotons, antideuterons, or neutrinos from 
dark matter annihilations in the galatic halo or the center of the Sun). A weak scale cross 
section can naturally lead to a dark matter density of the requisite value to explain the 
WMAP determination 6|]. Our interest here is in the tests of weakly interacting massive 
particles (WIMP) models via their elastic scattering on nuclei. 

In WIMP models a discrete symmetry exists that makes the DM candidate stable. In the 
minimal supersymmetric standard model (MSSM), the DM particle is the lightest neutralino, 
a spin 1/2 Majorana particle. In a non-SUSY context, a heavy Dirac fermion can be a viable 



DM candidate 
14, 3, 16, 17, 
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9|. Models such as Minimal Universal Extra Dimensions 
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191 ] incorporate a Kaluza-Klein (KK) parity which keeps the lightest KK 



particle stable. Here, the DM particle is spin-1. In the Littlest Higgs model solution to the 



gauge hierarchy problem, a discrete T-parity ensures that a spin-1 particle is stable 
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2l|. 



Scalar DM, with a spin-0 field can be a consequence of singlet Higgs field extensions o 



Standard Model which have been explored in a number of works 
and a model with two extra dimensions j^. There are a variety o: 
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particle that we do not specifically consider here, see e.g. 28 



29, 



other models with a DM 
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31 



32] 



SI scattering occurs though Higgs exchange and thus may well be much smaller than SD 
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scattering which occurs through Z-exchange 1 . The calculation of the SI cross section has 
substantial uncertainties from the estimate of the Higgs coupling to strange quarks 33l. |34|. 



The first WIMP searches focused on spin-independent (SI) scattering on a nucleus and have 
placed limits approaching 10 -8 pb. The SI experiments are designed to detect the coherent 
recoil of the nucleus caused by the DM scattering. For a heavy nuclear target, the coherent 
scattering increases the cross-section by the square of the Atomic Number. The best current 
limits on SI scattering have been placed by the CDMS 5 tower 35| experiment and the 



XENON10 experiment 36j. The CDMS upper bound on the WIMP nucleon cross section 
is 4.6 x 10~ 8 pb for a WIMP mass ~ 60 GeV. The XENON10 bound is 4.5 x 10" 8 pb for a 
WIMP mass of ~ 30 GeV and 8.8 x 10~ 8 pb for a WIMP mass of ~ 60 GeV. 

The present SI experiments have detector masses of order 10 kg. The next generation 
experiments will have about 100 kg mass. Detectors with 1000 Kg (IT) size are under 
design. The XENON100 experiment is expected to reach a SI cross section sensitivity of 
10" 9 pb for a 100 GeV WIMP Q- The LUX experiment, with a 100 kg xenon detector, 



expects to reach 4 x 10~ 10 pb at M DM ~ 40 GeV |38(. The ultimate goal of XENON1T 
is sensitivity to a SI cross section of 10~ 10 pb. One ton Liquid Argon detectors promise 



competitive sensitivities 39|, |40j, |41| 



The SD cross-section occurs through the axial vector coupling to the spin content of the 
nucleus; there is a J(J + 1) enhancement from the nuclear spin J. Until recently, the SD 
cross section limits had been about 6 orders of magnitude weaker than for SI, with the 
best bound from ZEPLIN-II at 0.07 pb for a DM particle of mass ~ 50 GeV scattering 
on neutrons 42j. The NAIAD, COUPP, and KIMS experiments had placed upper bounds 
on the cross section for SD scattering of dark matter on protons of 0.5 pb, 0.3 pb and 



0.2 pb, respectively, for a DM mass of order 100 GeV [43|, |44j, |45j. The SuperKamiokande 
(SK) search for neutrinos from DM annihilations in the Sun placed a stronger limit, albeit 
model dependent, on the SD cross section; converting the SK neutrino flux limit requires 
assumptions about the DM mass and the cross section for DM annihilation to neutrinos. For 
a DM mass of order 100 GeV the SK limits on the SI and SD cross sections are of order 10 -5 



and 0.6 x 10 2 pb, respectively 46|, |47|]. Liquid Xenon has about 50 percent odd isotopes 



1 In SUSY, sfcrmion exchanges can contribute to both SI and SD scattering. However, their contributions 
are typically small, even for light squarks. 
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which allows a measurement of SD scattering. New SD limits have been reported by the 
XENON10 experiment of 0.5 x 10 -2 pb for scattering on protons and 0.5 pb for scattering 
on neutrons at a WIMP mass of ~ 30 GeV 48j. Models generally predict the same cross 



section for scattering on protons and on neutrons, so the most restrictive of the two can be 
imposed in constraining models. 

New techniques with bubble technologies offer great promise for exploring much smaller 
SD cross sections in future experiments. The Picasso experiment in SNOLAB uses a su- 
perheated liquid in a gel. Its current proton- WIMP SD cross section bound is of order 1 
pb [49)]. In a future third stage of this experiment, with a 100 kg detector, sensitivity for 
SD scattering down to 10~ 4 pb is anticipated. The COUPP experiment 4J, [50] uses a 
continuously sensitive bubble chamber with compounds CF 3 I and C4-F10 of fluorine, F 19 , 
and iodine, J 127 . WIMP scattering from fluorine is determined by <jsd and scattering from 
iodine is determined by asi- 



y 



yield an observab 
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e signal of high 



561 ] . A large signal 



A strong SD cross section > 10~ 4 pb would also like . 
energy neutrinos from the Sun due to DM annihilations 
in neutrino telescopes such as IceCube of high energy neutrinos from DM annihilations in 
the Sun can occur only if there is a large SD scattering (since the SI cross section is already 
known to be <C 10~ 6 pb. The capture rate is \^ 



C Q = 3.4 x lO^s" 1 



Plocal 



0.3 GeV/cm 3 \ vi ocat 



270 km/s\ 3 /o-fo + crgj + 0.07<rf/\ / 100 GeV 



10" 6 pb 



m v o 

Al 



where pi oca i and vi oca i are the local density and velocity of relic dark matter, respectively. 
Present limits on the SI and SD scattering cross section on nuclei strongly suggest that an 
observation of these neutrinos from the Sun would only result from SD capture of WIMPS by 
hydrogen in the Sun. To calculate the expected signal rates, we closely follow the approach 



of Ref. 



56]. 



Recently, the DAMA/LIBRA experiment and the predecessor DAMA/Nal experiment 
have reported a combined 8.2a evidence for an annual modulation signature for DM particles 



in the galactic halo [58j. The implied SI cross section would be about 10 -6 pb. A possible 
way to reconcile the tension with other more restrictive bounds from other experiments is 
that DAMA is not detecting nuclear recoils, but electromagnetic energy deposition from 
axion-like interactions in the detector. For now, we set aside this interesting experimental 
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development and focus on nuclear recoil detection of DM. 

Our study is devoted to the analysis of what can be learned from the combined measure- 
ments of SI and SD processes, which can provide a powerful diagnostic in differentiating 
models. We specifically detail representative predictions of mSUGRA 2 , singlet extensions 
of the SM and the MSSM, a new Dirac neutrino, Littlest Higgs with T-parity, and Minimal 
Universal Extra Dimensions models (mUED). We assume that the DM particle is in thermal 
equilibrium in the early universe and we impose the requirement that the relic density of 
DM reproduces the measured WMAP density within ±10%. 

We also present general formulas for SI and SD cross sections for the scattering of WIMPS 
of intrinsic spin 0, 1/2, 1, and 3/2. 

The organization of the paper is as follows. In section [Til we outline the scan proceedure 
while in section [TTT1 we elucidate the DM scattering predictions of representative models that 
provide a WIMP candidate and detail how the models could be distinguished by measure- 
ment of SI and SD scattering. Then we proceed, in section IVIl to a model independent 
analysis of DM scattering for DM particles of arbitrary spin. Finally, we summarize the 
results of our study in section IVHI 



II. SCAN TECHNIQUE 

To arrive at the model predictions, we use a Markov Chain Monte Carlo (MCMC) ap- 
proach that is based on the Bayesian methods to scan over the relevant model parameters. 
With a MCMC, a random point, Xi in parameter space is chosen with a likelihood, Ci, 
assigned to it based on the associated model predictions such as the DM scattering cross 
sections and relic density. We then choose another random point, iCj+i, and compare its 
likelihood with the previous point in the chain. If the likelihoods satisfy 



> C, (2) 



where ( is chosen from a uniform distribution on the unit interval, the point is accepted 
and appended to the chain, otherwise a copy of x$ is appended to the chain. 
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Recently, efforts have been made to distinguish specific areas of mSUGRA parameter space that yield 



similar collider signatures using measurements from direct detection experiments [59 1 
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The power of the MCMC approach is that in the limit of large chain length, the distribu- 
tion of points, Xi approach the posterior distribution of model parameters for the constraining 
data. The speed at which the parameter space is scanned by MCMC is generally faster than 
traditional random or grid scans, with the difference in speed between the two methods 
increasing with larger parameter dimension. 

We construct the likelihood from observables which have a central value, fj,j, with uncer- 
tainty, (Tj, through the corresponding y 2 value 



(3) 



where dy is the model prediction for the j th observable of the constrained data based on the 
parameters in the chain Xj. To impose exclusion limits, we follow the approach of Ref. 60]. 

To construct the posterior distribution, we obtain the covariance matrix among the rel- 
evant parameters of the first 200 unique links in the MCMC chain and discard them as a 
"burn-in". We define the ranges over which the parameters are allowed to vary later in 
Section II III for each model. We use this distribution information to construct more efficient 
chain proposals for subsequent links in the new chain as outlined in Ref. 6l|. 



Convergence of the chain is tested using the method of Raferty and Lewis [62|, [63|, [64|. In 
checking that the chain has converged, we require that the posterior distribution is within 
2% of the 95% quantile with 95% C.L. This typically results in a few thousand unique links 
in the chain. We then double check convergence using the method of Ref. l65j. For more 



extensive reviews of this type of approach to parameter estimation, see Refs. [60|, 



61. 



65 



66|. 



III. PREDICTIONS OF MODELS FOR DM SCATTERING 

The measured DM density from the WMAP5 analysis of cosmological data (cosmic mi- 
crowave background radiation data combined with distance measurements from Type la 



supernovae and baryon acoustic oscillations in the distribution of galaxies) 



6| is 



Q DM h 2 = 0.1143 ±0.0034, h = 0.701 ±0.013 (4) 

We scan over the parameters of models with the constraint that the models give this relic 
density within ±10% assigned uncertainty associated with combined galactic and particle 
physics uncertainties. We separately consider model parameters that give a DM density 
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below the measured value to allow for the possibility that there may be more than one 
contributing DM particle. Our calculations of SI and SD cross sections and the relic density 



are made using the Micromegas 2.1 code |67j with an isothermal dark matter density profile 
and a relative DM velocity in the galactic halo of Vdm — 220 km/s. In calculating the SI 
scattering cross sections we adopt the recent determinations of the pion-nucleon sigma term 
a nN = 55 MeV and the sigma term that determines the mass shift of the nucleon due to 
chiral symmetry breaking <7o = 35 MeV 67J, [68 , 



69 



70]. 



A. mSUGRA model 

Supersymmetry (SUSY) has a conserved i?-parity that renders the lightest supersymmet- 
ric particle stable. Supersymmetry stabilizes the radiative corrections to the Higgs boson 
mass and realizes Grand Unification of the electroweak and strong couplings. Moreover, 
the quantum numbers of the Standard Model particles are explained by a singlet 16 rep- 
resentation of SO( 10) for each generation of fermions. The necessary breaking of the su- 
persymmetry is achieved in the minimal supergravity model (mSUGRA) by Planck scale 
mediation between the observed and hidden sectors of the theory. The mSUGRA model 

71 . 72I . I73I ]. It has only a small number of parameters 



has become a reference standard 
(m , mi/2, A , tan/3, sign(fi)) 3 . There are four regions of mSUGRA parameter space that 
give values of the relic density at or below the WMAP value. These are: 

• Focus Point Region (FP) [74], H, 76], also called the Hyperbolic Branch: This region 
is preferred by b — r unification and new experimental results from b — > 57 60]. 
The SUSY flavor changing neutral current and CP violating problems are resolved 
naturally by large sfermion masses. At large mo, the superpotential Higgsino mass 
term /i becomes quite small, and the lightest neutralino, Xi, is a mixed higgsino-bino 
state. Neutralino annihilation in the early universe to vector bosons is enhanced. This 
results in an enhanced SD cross section and enhanced neutralino annihilation in the 
Sun to neutrinos. The FP region would allow a precision gluino mass measurement 



The common scalar, gaugino masses and the soft trilinear terms are unified at the GUT scale to chosen 
values of m , m 1 / 2 and A , respectively. After specifying the ratio of the two Higgs VEVs, tan/3 = 
and the sign of the Higgsino mass parameter, /j,, the model at the electroweak scale is fully defined by the 
GUT scale parameters through the RGE running. 
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FIG. 1: Scan over the common scalar, m$, and gaugino mass, w-i/2) to satisfy relic density and 
LEP2 constraints on mSUGRA with specific values of Aq and tan (3 (A$ = 0, tan f3 = 30, 55). Open 
squares (in red) show parameter values that have a relic density within ±10% of the measured value; 
solid points (in black) have lower relic density values. 



(±8%) for Mg in the range of 700 to 1300 GeV 77j] - The FP region falls into a more 



general class of neutralino models which are termed "Well-tempered Neutralinos" in 
which the lightest neutralino is a Bino-Wino or Bino-Higgsino mixture H]- 

A- Funnel Region (AF): The A-funnel region occurs at large values of the parameter 
tan/3 pa 50, near 2M x o ~ m^. The neutralinos annihilate through the broad pseu- 
doscalar Higgs resonance A 80, 81, 8^|. There is also a light Higgs resonance region 



where 2M x o ~ m^ at low values 



S3 



84]. 



Coannihilation Regions (CA): The neutralino-stau co-annihilation region occurs at 
very low mo but any jn-i /o values, so that ~ M x o, and neutralinos can annihilate 
against tau sleptons 85j, |86| in the early universe. For certain A Q values which dial 



mi to very low values, there also exists a stop- neutralino co-annihilation region 871 ] . 



Bulk Region (BR): The bulk region is at low rriQ and low m 



where neutralino 



annihilation is enhanced by light t-channel slepton exchange |81l . [84J. The WMAP 
determination of the relic density has pushed this allowed region to very small mo and 
mi/2 values, while LEP2 limits on M ± and m^ exclude these same low values so that 
the bulk region is disfavored [88J, |8jj, [90]. 



8 



mSUGRA 




FIG. 2: Neutralino relic density and LEP2 allowed regions in the mSUGRA model obtained from 
a scan over the common scalar, mo, and gaugino mass, min, while allowing variations of the 
parameters Ao,t&nf3 and mt- One, two and three sigma contours are shown 



We note that neutrino Yukawa couplings can significantly affect the RGE evolution and 
the calculated neutralino relic density in regions of parameter space where soft SUSY- 
breaking slepton masses and / or trilinear couplings are large 9l| . The changes can be large in 
the focus point, A-funnel, and stop-coannihilation regions of mSUGRA and can expand the 
allowed regions of the mSUGRA parameter space and the dark matter direct and indirect 
detection rates. We do not fold in the uncertainties associated with the neutrino couplings 
in the present study. 

Figure [1] shows representative regions in mSUGRA parameter space where the relic den- 
sity is accounted for within ±10% of the WMAP value (red open points) or is less than the 
WMAP value (black solid points). The choices of Aq = 0, fi > 0, m t = 170.9 GeV and 
tan j3 = 30 (Fig la) and tan j3 = 55 (Fig lb) are chosen to illustrate the FP, CA and AF 
regions. Well delineated regions that satisfy the relic density are carved out of the param- 
eter space. As noted above, these regions become wider when the unknown masses of the 
right-hand neutrinos are taken into account in the evolution from the GUT scale 91]. 

When we allow the trilinear parameter, A , and tan/5 to vary, as well as the input value 
of the top quark mass mt, we arrive at a broader range of parameter values that reproduce 
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the DM density. In our analysis, we allow the MCMC to scan within the parameter ranges 

50 GeV < m < 5 TeV 

150 GeV < m 1/2 < 1.5 TeV 

1 < tan/5 < 50 

~ " (5) 

-3 TeV < A < 3 TeV 

< fi 



We allow a gaussian variation of the top quark mass around its central value of m t = 170.9 
GeV with a standard deviation of 1.8 GeV 92|. The marginalization over the "nuisance 



parameters" allows the uncertainty of the input parameters to be reflected in the posterior 
distributions of the model parameters. In mSUGRA, the top quark mass variation strongly 
affects the running from the GUT scale to the weak scale. Therefore, through its variation, 
the top mass uncertainty expands the parameter ranges that are consistent with the observed 
relic density compared to results obtained by fixing the top mass as is often done in such 
studies. The resulting distribution of the common scalar and gaugino masses (m ,mi/ 2 ) is 
shown in Fig. El 

In general the neutralino composition in terms of the gaugino (B, W 3 ) and Higgsino 
(Hi, H 2 ) states is given by 

Xi = N n B + N U W 3 + 7Vi3#i + N 1A H 2 (6) 

The spin dependent scattering cross section is largely governed by Z-boson exchange and is 
sensitive to the Higgsino asymmetry 

<t 3 d<x (7) 

The FP region is unique among the mSUGRA regions that reproduce the relic density in 
having a large Higgsino asymmetry and a correspondingly large a$D 

The SD and SI scattering cross sections in mSUGRA are displayed in Fig. [31 where it 
is apparent that different solutions to the DM relic density populate different regions of 
<Jsd versus a si- The limit of current DM recoil experiments is strongest in the mass range 
Mdm ~ 50 — 100 GeV. In all subsequent a si vs. ctsd results, we present the current or 
projected best limit of the recoil experiments. The use of <Jsd/&si to distinguish mSUGRA 
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FIG. 3: Posterior distributions of SI and SD cross sections in the mSUGRA model. The region 
of high SI and SD cross sections corresponds to the FP region and should be probed fully by the 
XENON10 and Super CDMS 25 Kg experiments. The tail that extends to lower SI and SD values 
corresponds to the AF and CA regions. The proposed COUPP1T experiment should probe these 



regions 



regions has also been advocated previously, see e.g. Ref. [93J. However, it should be noted 
that the experimental limits are progressively less constraining at higher and lower M^m 
than the best limits. The FP region with its relatively large <jsd can be definitively tested 
by SD, SI measurements. DM detection corresponding to the FP region would have major 
significance for colliders in that high mass sfermions would be implied. Similar to the spin-0 
DM case, detection of SD scattering of DM would immediately rule out a model with a 
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gaugino dominated neutralino 



94J which predicts a vanishingly small SD cross section. 




FIG. 4: Posterior distributions of the induced muon flux at the Super-Kamiokande detector 
(cm _2 s~ 1 units) from high energy neutrinos created by DM annihilations in the Sun for mSUGRA. 
Only the low neutralino mass FP region is constrained by the Super-K limit. 



The Su per -Kamiokande experimental upper bound on the induced muon flux is 5 x 10 



— 2 — 1 

cm s 



-15 



47], and is shown by the horizontal line in Fig. HI This bound only eliminates 



only a small fraction of the FP parameter space. 



B. Singlet Higgs Extensions of the MSSM 



The motivation for appending a singlet Higgs field S to the MSSM is to explain why the 
Higgsino mixing parameter \x in the superpotential 



W = nH u - H d 



(8) 



is at the TeV scale. The dynamical \x solution is to let the \i parameter be generated by a 
vacuum expectation value of a S of order Msusy- 



/icff = X(S) 

Th ree c l asse s of s inglet model e xten sions xMSSM have been considered 



(9) 
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100 
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102, 
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95 



96 
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104 



105 



98, 



106J with additional Higgs-sector superpotential terms as 



follows: 



12 



(i) NMSSM: S 3 (cubic) 

(ii) nMSSM: S (tadpole) 

(iii) UMSSM: U{1)' symmetry 

In the xMSSM the spin-0 S field mixes with the neutral Higgs fields of the MSSM; the spin- 
1/2 singlino field S in the composition of the lightest neutralino mixes with the neutralinos 
of the MSSM. The consequen ce for DM is that the S component changes the DM couplings 



and cross sections | 100 
from the MSSM [10 



104l | . Th e predictions of the tadpole model are the most different 



104 Il05| : we therefore focus on it here. The singlet scalar field 



also affects LHC phenomenology as the S admixture in the Higgs states changes the Higgs 



production and decay r ates and modifies the cascade decay chains of the sparticles 



98 



100 



101 



105 



107 



108]. 



99, 



We scan over the parameters directly relevant to the DM phenomenology of this model. 
Their ranges are chosen to be 



1 


< 


tan (5 


< 


20 


200 GeV 


< 


(S) 


< 


2 TeV 


100 GeV 


< 


Veff 


< 


1 TeV 


50 GeV 


< 


Mi 


< 


1 TeV 


100 GeV 


< 


M 2 


< 


1 TeV 


100 GeV 


< 


A x 


< 


1 TeV 


100 GeV 


< 


A x 


< 


1 TeV 


(400 GeV) 2 


< 


t F 


< 


(400 GeV) 2 


(200 GeV) 3 


< 


t s 


< 


(100 GeV) 3 



(10) 



where the parameters A\,tp and ts are defined in Ref. 1001 ] . We also place a minimum cut 
on the lightest neutralino mass of 10 GeV since the current DM scattering experiments have 
considerably reduced sensitivity to a DM mass this low. 

Figure [5] shows the SI and SD cross sections for the nMSSM. The contours show the 
posterior distributions with the present SI and SD experimental constraints enforced. The 
SD cross section can be rather large since the dominant annihilation occurs through the Z 
boson. Therefore, to counter the small annihilation rate in the early universe due to the 
small neutralino mass, the neutralino pair must have a larger Z boson coupling. 
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FIG. 5: Posterior distributions of SI and SD cross sections that satisfy the relic density and LEP2 
constraints in the tadpole extended SUSY model. This model has light neutralinos with relatively 
large SI and SD cross sections that make the prospects for discovery promising. 



The nMSSM predictions of a si and a$D will be largely probed by future experiments, 
provided that the neutralino mass is not too small. Relaxation of the 10 GeV lower cutoff 
on the DM mass would be to loosen the prospects of DM detection in this model. 



C. Scalar singlet extended SM (xSM) 



One of the simplest models that includes a viable DM candidate is a real singlet Higgs 



extended SM with a discrete Z2 symmetry 23|, |24j, |25|]. In this model, a stable singlet is 



appended to the SM and interacts with the SM fields only through the Higgs boson. As 
a consequence, the relic density and elastic scattering predictions are correlated with the 
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Higgs sector. The potential for a stable singlet Higgs boson, S, is given by 



(11) 



where H is the SU(2) doublet field and m 2 , A are the usual SM parameters of the Higgs 
potential. Combining this singlet extended Higgs sector with the rest of the SM gives the 
"xSM", the extended Standard Model. The singlet DM mass is determined by M^ M = 
\k 2 + \5 2 v 2 , where v = 246 GeV is the SM Higgs VEV. In the MCMC scan, we allow the 
parameters to vary within the ranges 



114 GeV < M h < 150 GeV 

50 GeV < M DM < 1 TeV 

< 6 2 < 4 

< k 4 < 4 



(12) 



where the upper limit on the Higgs mass comes from electroweak precision measure- 
ments 251 ] . 



Further extensions of this model include making the singlet a complex field 109| or adding 
right handed n eutri nos, which would explain neutrino mass generation when the singlet 
obtains a VEV llOj. The real singlet Higgs field has also been shown to enhance the first 



order phase transition of electroweak symmetry breaking that is necessary for electrweak 
baryogenesis and may reduce the tension between electroweak precision observables 111 ] 
and the non-observation of the SM Higgs boson below 114 GeV 25]. 
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FIG. 6: Real singlet extended SM posterior distributions of the SI cross section vs. Moai- 
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Since the DM candidate in the xSM class of models is spin-0, the SD cross section 
vanishes. If a SD signal is seen, this class of models would be ruled out. The SI cross section 
is generally small, below ~ 10~ 8 pb as seen in Fig. [61 

The xSM scenario could be indicated by an observation of a SI cross section of O(10~ 9 pb) 
with next generation detectors and a null result in searches for SD scattering and high energy 
neutrinos from the Sun in neutrino telescopes. 



D. Dirac neutrino DM 



A simple DM model can be realized with a stable neutral heavy Dirac fermion 112 . 



113j |. The stability of the heavy neutrino can be guaranteed by imposing baryon number 
conservation in a warped GUT model. In this model, the heavy neutrino can be 0(100 GeV) 
and annihilate in the early universe via the Z, Z' and h bosons. The parameter ranges over 
which we allow the MCMC to scan are 



114 GeV 


< 


M h 


< 


1 TeV 


10 GeV 


< 


M v , 


< 


1 TeV 


300 GeV 


< 


M w i/ Z > 


< 


2 TeV 





< 


9h 


< 


0.5 





< 


9z 


< 


0.5 





< 


gz' 


< 


0.5. 



(13) 



where Mh,M u i and Mw'/z' are the Higgs boson, Dirac neutrino and W'/Z 1 gauge boson 
masses, respectively. The couplings of v 1 to the Higgs, Z and Z' bosons are gh,gz and gz>, 
respectively. Since the Higgs boson couples the neutrino and matter with Yukawa strength, 
it often has an inconsequential effect on the relic density and scattering rates. The Z and Z' 
gauge bosons have gauge strength interactions, and the SM Z boson has the strongest effect 
due to its lower mass. In this scenario, where the Z boson dominates in the calculation of 
both the relic density and elastic scattering cross section, the SI and SD cross sections are 
tightly correlated, as seen in Fig. [7J 

This model predicts a rather large SI cross section that is partially excluded by limits from 
XENON10 and CDMS 2008. The contours in Fig. [TJshow the posterior distributions with 
the present SI and SD constraints enforced. The predicted neutrino flux from annihilations in 
the Sun is significantly smaller than if the scattering limits were not included. Nevertheless, 
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FIG. 7: Posterior distributions of SI and SD cross sections in the Dirac neutrino model of Ref. 
This model has a heavy stable neutrino that annihilates predominantly through Z and Z 1 gauge 
bosons. 



the prospects for a significant signal in IceCube from this DM candidate are still good, as 
detailed in Section [V] 



E. Minimal Universal Extra Dimensions 



12 



xtra d imension models have gained popularity as a viable beyond the SM alternative 11 
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ia 



114 



1151 ] . Many types of extra dimensional models have been proposed. Here we 



focus on the minimal Universal Extra Dimensions model (mUED) 121 ] . The stability of the 
DM candidate is ensured by KK-parity, which results from momentum conservation in the 
extra dimension. 

In the mUED model the DM has spin-1; it is the KK partner of hypercharge gauge boson, 
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B\. The predictions of the KK particle spectra are mainly determined by one significant 
parameter, the curvature of the extra dimension, R. Additionally, radiative corrections from 
the higher KK levels shift the masses of the KK states in the first level. The number of KK 
levels included in the mass renormalization is given by AR, where A is the cutoff scale of the 
model. Many KK states from the first excitation in this model are nearly degenerate with Bi, 
requiring inclusion of coannihilation processes in the calculation of DM relic density. After 
the relic density constraint is imposed, the curvature is found to be about R" 1 ~ 600 — 700 
GeV. We allow the MCMC to vary over parameter ranges 

114 GeV < M h < 600 GeV 

400 GeV < R- 1 < 800 GeV (14) 
5 < AR < 100. 



where the upper bound on the Higgs mass 19j, Ill6j is from electroweak precision constraints. 



The DM scattering cross sections in mUED are displayed in Fig. [HJ This model has a 
predicted 'sweet spot' of asD ~ 0(1O~ 6 ) pb for which the relic density is reproduced; the 
<Jsi range is considerably dispersed in comparison. The source of this narrow range for the 
SD cross section is the limit from the relic density on the curvature of the extra dimension. 
Since the KK quarks have approximately the same mass as the inverse curvature, the SD 
cross section is closely tied to the relic density. The SI cross section is more dispersed due 
to the relatively large variation of the Higgs boson mass. The test of the direct detection 
predictions of this model may require a 1 Ton detector. However, an observation of a 
neutrino flux resulting form B\ annihilations in the Sun is possible: see Section [V] 



F. Littlest Higgs with T-parity 



In Little Higgs models, the Higgs is an a ppro x imate Go l dsto ne boson with its mass 



117 



118, 



119 



120|. Additional states are 



protected by approximate global symmetries 
introduced to preserve this global symmetry and keep the Higgs from receiving quadratically 
divergent radiative corrections. The Littlest Higgs model is one of the simplest and popular 
models 20j, |2lJ. It has a SU(5)/SO(5) symmetry breaking pattern which gives two gauged 
copies of SU{2) x U{1) that are subsequently broken to SU{2) L x U(l)y at the scale / ~ 1 
TeV. 
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FIG. 8: mUED predictions of the SI and SD cross sections. The mUED predictions for the SD 
cross section may be probed with a one Ton detector if it has sufficient sensitivity at Mom ~ 700 
GeV. 



While the Littlest Higgs model protects the Higgs boson mass, it does not satisfy elec- 
troweak precision constraints. Large corrections to electroweak observables occur at tree- 
level. This may be allev iated by i mpo sing a discrete T-parity that interchanges the two 



copies of SU(2) x U(l) [121 



122 



1231 ] . Once T-parity is introduced (the LHT model), 



electroweak corrections scale logarithmi cally with the cutoff scale, A ~ Airf. We assume 



that the Wess- Zumino-Witten anomaly 
candidate 1251 ] . 



1241 ] does not break T-parity to have a viable DM 
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The parameter ranges that we allow for the LHT model are 

114 GeV < M h < 1 TeV, 
500 GeV < / < 2 TeV, 



(15) 

2 < tana < 2, 
< k, hi < 4, 



where the k, ki couplings in the Lagrangian (see Eq. 2.14 of Ref. [1261 ]) determine the masses 
of the T-odd quarks and leptons masses, respectively. The angle a defines the mixing 
between the T-even heavy quark and the SM top quark. We require the breaking scale, 
/, to be below 2 TeV to avoid a reintroduction of the hierarchy problem. In additio n, we 



included the electroweak precision constraints on the oblique parameters given in Ref. 126] 



and calculated the additional \ 2 contributions associated with these parameters according 



to Ref. 25[. The electroweak constraints in this model can allow a larger Higgs mass than 



in the SM, which will result in a suppression of the SI cross section. To satisfy experim ental 



bounds on four-fermi operators, the values of the T-odd masses must satisfy the bound 

M^odd 



126] 



GeV 



< 4.8 x 10- J (//GeV) 2 (16) 



Associated with the introduction of T-parity, there is a stable spin-1 DM candidate. All 
non-SM gauge bosons and the triplet Higgs are T-odd, while SM fields are T-even. The 
lightest T-odd state is the heavy photon, Ah- The relic density is quite sensitive to the 
Higgs mass as the heavy photon annihilates dominantly through the Higgs boson. Other 
annihilation processes involving the new heavy quarks are generally suppressed their large 
masses and the small hypercharge, Y = 1/10 of the DM particle. The SI scattering proceeds 
through the Higgs and T-odd quarks, giving SI cross sections of order O(10~ 9 — 10~ 10 ) pb. 

We present the SI and SD cross sections in Fig. [91 Since SD scattering of the DM heavy 
photon occurs only via the T-odd quarks, the SD cross section is suppressed. However, the 
predicted range of the SI cross section is within the projected sensitivities of future one Ton 
experiments. A modification of the LHT model has bee n re cently proposed that has the 



1271 ]. The coupling enhancement 



DM candidate couple more strongly to the Higgs boson 
is about a factor of 10 larger than that of the standard LHT model and the model yields a 
much higher SI cross section. 
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FIG. 9: Littlest Higgs model with T-parity predictions of the SI and SD cross sections. The SD 
scattering is highly suppressed. 



IV. DM SCATTERING UNCERTAINTIES 



Significant uncertainties on the hadronic matrix elements [33J are associated with the 
values of a n N an d <Jo'- c.f. Section HLT1 The dominant contribution to SI scattering is Higgs 
exchange coupled to the strange quark. The uncertainty on the strange quark contribution 
causes variations in the SI cross section by an average of about ±30%. In the left panel 
of Fig. dni we present the posterior distributions of the SI and SD cross sections in the 
mSUGRA model with the uncertainties in the sigma terms neglected. In the right panel of 
Fig. [10] we include the uncertainties on the hadronic matrix elements. These uncertainties 
substantially expand the regions in the SI-SD plane encompassed by the predictions. In all 
other illustrations in this paper we include the hadronic uncertainties. 
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mSUGRA: SI vs. SD, Hadronic uncertainties not included 
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FIG. 10: Uncertainties in the sigma terms can produce substantial uncertainties in the SI and SD 
cross section predictions. We illustrate this effect in the mSUGRA model; the left panel neglects 
these uncertainties and the right panel includes the uncertainties. 



V. EXPECTED SIGNAL RATES AT ICECUBE 



High energy neutrinos are expected from t 



gravitational capture in the Sun 51 



3 



ie annihilations of DM that accumulated via 



128, 



1291 ] . Neutrino telescopes are well positioned 



to search for neutrinos of this origin 53]. The IceCube experiment at the South Pole is 
underway and expects to have about 50,000 events from atmospheric neutrinos in the near 
future [5J, 1 130t ] , demonstrating the capabilities of the detector to identify neutrino events. 

The basic requirements for a strong signal in km 2 sized neutrino telescopes from DM 
annihilations in the Sun are as follows: 

• A relatively large DM annihilation rate either directly to neutrinos or to SM parti- 
cles that subsequently decay to high energy neutrinos (gauge bosons, top and bottom 
quarks, and tau leptons). The latter is generically the case for all the models consid- 
ered. 

• The DM mass exceeds the muon energy threshold of the neutrino telescope, which is 

~ 50 GeV for IceCube 4 . Since the DM annihilation proceeds non-relativistically, 
the DM mass is the maximum energy that any neutrino can have. This kinematics 



A future inner ring of photomultiplier strings in IceCube may lower this energy threshold to 10 GeV. 
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IceCube muon rate 



IceCube signal significance 




FIG. 11: The expected numbers of events and the statistical significance of DM signals for 3 and 5 
years running of the IceCube neutrino telescope. Only the FP region in mSUGRA and a portion 
of the mUED parameter space give a significance > 5er for 5 years of running. 



essentially eliminates current IceCube tests of the tadpole nMSSM model since the 
lightest neutralino mass is typically only 20-30 GeV, although in some small, fine- 
tuned portions of parameter space, the DM mass may extend above the Z-resonance. 

• A sufficiently large SD cross section for neutralino captures (see Eq. [I]) to produce a 
detectable neutrino flux. A SD cross section > 10~ 4 pb is necessary to yield a neutrino 
flux that could be observed by IceCube. Models that have either no or a vanishingly 
small SD cross section have no chance of producing an observable signal in km 2 sized 
neutrino telescopes. This includes the xSM and the Littlest Higgs with T-parity. 

The predictions for IceCube of the models with SD scattering are given in Fig. [II] The 
numbers of muon events and the corresponding statistical significances for the models are 
given for both 3 and 5 years exposure with the full array (50% of the array is now in place). 
The calculations follow Ref. 56|. The FP region of mSUGRA bodes well for the detection of 
these neutrino events; the predictions of the other models are below the 5a discovery level; 
the RHN and mUED models are on the borderline of 3a detectability. 



VI. MODEL INDEPENDENT DM SCATTERING CROSS SECTIONS 

In this section we categorize the interactions of DM candidates according to their intrinsic 
spins and provide generic expressions for the SI and SD cross sections. 
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A. Self- Conjugate Spin-0 



A self-conjugate Spin-0 DM particle, 0, cannot have a SD scattering amplitude at thresh- 
old. The generic effective Lagrangian is 

C* 3 -£4>4> q(S, + iP, 7 5 )9 ■ (17) 

The interaction, flipping chirality, is proportional to the q nar k mass rn,. Other forms of 
the interaction can be reduced to this form by the equation of motion [27[ with an emerging 
mass m q factor. Since the scalar DM candidate is self-conjugate, the vectorial part of the 
Lagrangian vanishes 5 . The pseudoscalar coupling, which is derivative in nature, is not a 
leading contribution to the scattering because the scattering amplitude is proportional to 
the DM velocity. The surviving scalar interaction S is purely SI. For a nucleus target (N) 
of Z protons and A — Z neutrons, the accumulated amplitude is proportional to the nucleus 
factor Cn, 

C N = ^(S c + S b + S t )Af TG + S ^ Z fr q + ( A ~ z )fr q ) • (18) 

q=u,d,s 

Here the first grouping is the gluonic contributions induced by the heavy quarks c, b, t, and 
the second grouping is the contributions of the light quarks, where is the fraction of 
the proton mass that the light quarks (u, d, s) represent, fT q m p(p\p) — ip\ m gQQ\p) with 
fx g = 1 — as the remaining gluon fraction. Their values for the proton and the neutron 



are 



131] 



fP u = 0.020 ± 0.004 , f% = 0.026 ±0.005, /£, = 0.118 ± 0.062, (19) 
/P u = 0.014 ± 0.003 , f% A = 0.036 ± 0.008 (20) 

The cross section for the scattering of a scalar dark matter particle on a nucleus target is 

"47 = ^-J^r-) 2 (^cS . (21) 



47rA 4 ym^ + mjvy 

The formulas above (and also below) for the cross sections are understood to be valid in the 
low momentum tranfer limit. A form factor can be included for the finite size of nucleus. 
For a proton target = m p , and A = Z = 1. 



This is in close analogy that the self-conjugate Majorana fcrmion cannot have a vectorial coupling. 
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If DM is spinless, there is no amplitude for SD scattering and no amplitude for static DM 
pair annihilation into a light fermion pair. For this reason, this model predicts no mono- 
energetic neutrinos 



55] from annihilations of DM in the solar core. However secondary 
neutrinos with energy distributions arising from DM to WW, ZZ, tt etc. are still allowed, 
as studied in [56]. 



B. Dirac Fermion 

If the DM particle is a Dirac fermion, the generic interaction is described by 
vectorial/axial- vector and scalar/pseuodoscalar couplings, 

£</, D — ^(S^, + rfP^ip qm q (S q + P q irf)q 
+^ 4>1»{V + A 7 5 )^ q-f(v q + a q ^)q 

The vector- vector coupling Vv q and the scalar-scalar coupling S^S q give rise to the spin- 
independent (SI) low-energy scattering amplitude. For a Dirac dark matter particle ip, the 
SI cross section for scattering on a nucleus N is 

d"s? = (^) (U) \ Vv u(Z + A) + Vv d (2A -Z) + S^m p /A)C N \ 2 • 

In this formula we have used valence quark model relations. Note that if the incident DM is 
an antiparticle ip, the cross section is given by the same expression except that the V term 
flips sign. Thus, particle i/j and the antiparticle ip would scatter on a nucleus with different 
cross sections when both vector- vector and scalar-scalar amplitudes exist. 

As SI measurements in ip + N scattering highly constrain the product Vv q , we shall 
neglect v q for simplicity and focus our attention on the possibly much larger SD scattering. 
The vectorial part V still plays a dominant role in the DM annihilation but not in the SD 
ifj + N scattering because its time-component does not match the space-component of the 
quark axial current. The SD differential cross section can be written as 



4 . , . ml 



da^ = -J(J+1 



7T A 4 



2 

dVt 2 / m^rriN 



™%n = 1 • (22) 

47r v \m^ + m Ni 



where J is the angular momentum of the target nucleus. In case of the proton target that 
is relevant for solar capture of the DM, J = |. The fractional quark-spin coefficient X g is 
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defined in the Appendix. All participating flavors q are summed at the amplitude level. The 
reduced mass of i/j + N appears as m^N. A direct search for DM is optimal with a target 
nucleus of comparable mass to the DM particle and with a large J value. A multiplicative 
unwritten form factor F(Q 2 ) due to the nucleus size describes the dependence on the squared 
momentum transfer, Q 2 . The Q 2 dependence may be significant for a large reduced mass. 

If the DM couples to the Z boson, the interaction can be written as C D ip'-f f j i (V+A'-f 5 )ijjZ fl . 
This effective amplitude can arise from the t-channel exchange of Z with the following 
relation 

a q _ 9 rp V q 9 frp ■ 2 Q \ 

A 2 2 cos OwMg A 2 2 cos 6wMz 
where T 3 u = ±, and T 3 d = -\ 



Static annihilation of ipip to a light fermion pair can take place through the vectorial 
amplitude V. The vv final state is of interest as a source of mono-energetic neutrinos jsol] 
from the core of the Sun. From DM annihilation through a Z boson, or a new gauge boson 
of a larger gauge group under which DM is charged, the fraction of the DM annihilations 
to neutrinos can, in principle, be quite large. Models that allow th e possibility of Dirac DM 



include the hidden fermion in the Stueckelberg Z' model [112|, [132|. Heavy sterile neutrinos 
also belong this type of model; they constitute warm dark matt er and are not considered in 
our study as we assume the DM velocity is nonrelativistic |l33 |. 



C. Majorana spin— | Fermion 

For a spin | dark matter (DM) Majorana fermion, the relevant dynamics is determined 
by the interaction 

£ D ^ x{S x + H 5 P X )X qm q (S q + P q i^)q + ^(x^TsX^K + a qlb )q . (23) 

There is no vector- vector interaction for the self-conjugate x field. The SI cross-section is 
given by the scalar-scalar interaction. The asi formula is like what is found in the Dirac 
case, with V = and an overall factor of (2!) 2 due to two counts of Wick contraction for 
the Majorana field, 
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The SD scattering cross section is given by 



dcr x+N 

aa SD 



(2<) 2 ^ + l)^ 



dQ 

47T 



(25) 



This result agrees with Eq. (2.26) of Ref. 1341 ]. provided that one is careful in distinguishing 
the amplitude and the Wilson coefficient. 

The annihilation of a pair of Majorana fermions \X & t res t a neutrino pair is forbidden 
because of the helicity suppression from the extremely small mass of the neutrino. The same 
helicity suppression applies to any light fermion pair. The helicity suppressio n ma y be lifted, 
however, if a Z boson is radiated from an internal or external neutrino line 1351 ] . However, 
this is a very weak process 0{b\a 2 rn^l A 4 ) and it would also distort the neutrino line shape 
consi dera bly. Models that allow this possi bility include SUSY {71, 72, 73], singlet femion 



DM 



1361 ]. and heavy Majorana neutrinos 1371 ] . For a detailed study of neutrino signals 



from neutralino annihilation in SUSY through intermediate states such as W and Z bosons 
and top quarks, see Ref. [56]. 



D. Spin-1 



13 



A spin-1 DM 



14,05 



3,Jl8 



particle is encountered in Universal 
, Littlest Higgs with T-parity (LHT) 



xtr a Dimens ions (UED) 



118 



120 



[lq, 



11 



12 



1381 ]. and Randall- Sundrum 



(RS) models 11], with an effective lagrangian 



C D S q qm q q + ^ ie ^{B u id a B,) q lpl5 q 



(26) 



The SI cross section has the same form as that of scalar DM 

(2!) 2 



do™ 



2 / \ 2 

w 



m B 



'N 



(27) 



47rA 4 + rug t 

The SD scattering amplitude B(ei) + N — > B(ef) + N is given in terms of the incoming 
(outgoing) polarizations e, (e/) by the expression 

M = 2\^(e}),te°^(e l ) u (2m B ) (N\q lpl5 q\N) 2m N (28) 

where the factor 2m N links the nonrelativistic convention {N\q^ b q \N) = 2\ q N) (N\J N \N) 
to the relativistic normalization. The combinatorial factor 2! is due to the self-conjugate B\ 
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field in the Lagrangian operator. The amplitude simplifies to 



M = ^A q (2m B ) 2m w (e f );^(6 i ) J • (N\q lkl5 q\N) 
The spin-1 matrix is (Sk)tj = itekj- Therefore, the amplitude takes the form 

M = ^(2m B ) (2m N ) A q 2X^S ■ J . 



A 2 



(29) 



(30) 



The SD scattering cross section is given by 



mi 



(2!)2 __ J(J+1) ^ 



dQ 



(31) 



Like fermionic DM, the SD operator does not give a static spin-1 DM annihilation amplitude 
to a fermion pair. However, there are other effective operators, as outlined in the Appendix, 
that are relevant to fermion final states. Mono-energetic neutrinos can be produced from 
spin-1 DM annihilation in the solar core. 



E. Spin | 



For the DM particle that is described by a self-conjugate Rarita-Schwinger field of 
spin S = §, like the gravitino G° in supergravity (SUGRA), the relevant interaction can be 
written as 

C D ^,S^ qS q m q q + ^l^tf^^fa^q) . (32) 
The SI formulas are like those of the spin-i Majorana case. 



da 



SI 



(2!) 



2 

2 m ^N 
VTA 4 



dVt\ 
4tt) 



m r . 



S^—^-Cn 



The SD amplitude in our special case S = 3/2 turns out to have a universal form, 
M(V + N -> * + N) = jA q (2m*) {2m N )S ■ {2\ q N) J) . 



(33) 



The square is simplified with the use of the identity Trf^S" 6 ] = 5 ab (2S + 1)S(S+ l)/3 when 
the spin orientation of S is summed. The spin-averaged square of the elastic SD scattering 
amplitude is 

2 



z — 'm s 



fM 2 m 2 N J(J + 1)(2S + 1)S(S + 1) 



(34) 
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Averaging over the (2S + 1) spin components, we obtain 



da SD (HN) = ^S(S + 1) J(J + 1)^ 



2 



^ IW 2 )| 2 (35) 



With appropriate substitutions, this generic formula agrees with that of the self-conjugate 
spin-1 case, as well as that of the Majorana spin-| fermion, if we set Aa q = \A q . 

VII. CONCLUSIONS 

Both direct detection experiments for DM scattering and LHC experiments can discover 
a weakly interacting DM particle and reveal its properties. In Section III II we showed that 
the SD versus SI cross section provides an excellent DM diagnostic. We highlighted this by 
scanning over the parameter space of six Beyond the SM models which supply viable DM 
candidates: mSUGRA, singlet extended SM and MSSM, a stable Dirac neutrino, Littlest 
Higgs with T-parity and Minimal Universal Extra Dimensions. To thoroughly scan the pa- 
rameter spaces, we adopted the Bayesian method that is the foundation of the Markov Chain 
Monte Carlo approach. The DM models can have distinct phenomenological predictions. 
We showed that the DM model possibilities can be narrowed by measurements of both SI 
and SD elastic scattering. The direct signals for DM in recoil and neutrino telescope exper- 
iments are complementa ry to LHC experiments in distinguishing the beyond the standard 



model physics scenarios 
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We summarize below the model predictions for the DM cross sections; the posterior 
distributions are summarized in Fig. [T2j 

• In mSUGRA, the FP region provides the largest SI and SD cross sections. This is due 
to the mixed Higgsino nature of the lightest neutralino; the neutralino couplings to 
the Higgs and Z bosons are large. The Bino nature of the lightest neutralino in the 
CA and AF regions causes these scenarios to have substantially smaller cross sections. 

• The tadpole nMSSM model has large SD scattering, of order 10~ 3 pb, and a wide range 
of SI cross section. This is a consequence of the DM annihilation occuring through 
the Z boson . To counter the small annihilation rate in the early universe (due to the 
small neutralino mass in the model), the neutralino pair is required to have a larger 
Z boson coupling, resulting in a large SD cross section. 
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FIG. 12: The posterior distributions of c^ A/_p and &^[f~ p for six models. 

In the singlet extended SM, the DM candidate is spin-0, which gives a vanishing SD 
cross section. The SI cross section is generally small, below ~ 10 -8 pb, and occurs 
through Higgs boson exchange. If SD scattering of DM is observed, the class of models 
with spin-0 DM would be immediately excluded as being the sole origin of the DM in 
the universe. 

For stable Dirac neutrino DM, the Z boson dominates in the calculation of both the 
relic density and elastic scattering cross section and makes the SI and SD cross sections 
tightly correlated and large. 

In mUED, a 'sweet spot' of <j,sd ~ O(10~ 6 ) pb exists for which the DM relic density 
is reproduced. The relic density is strongly dependent on the curvature parameter 
and fixes its value. The KK quarks have approximately the same mass as the inverse 
curvature and the SD cross section is thus closely tied to the relic density. On the 
other hand, the asi cross section is more dispersed due to the larger variation of the 
Higgs boson mass. 

In the LHT model, the SD interaction occurs through T-odd quarks which have a 
small hypercharge. Therefore, the SD cross section in this model is typically very 
small. In contrast, the SI scattering proceeds through the Higgs and T-odd quarks, 

30 



giving experimentally accessible SI cross section values. 




mSUGRA nMSSM xSM RHN mUED LHT mSUGRA nMSSM RHN mUED 



FIG. 13: Box and whisker plot of the SI and SD cross sections. The box represents the coverage 
of the middle 50-percentile and the whisker is ~ 100% coverage. The xSM is not shown in the SD 
panel as its cross section vanishes. The large SD cross section of the nMSSM is difficult to probe 
because of the low (~ 20 — 30 GeV) mass of the DM particle in this model. 

We provide a summary of the SI and SD cross sections by the box and whisker plots 
in Fig. [T3l The boxes represents the coverage of the middle 50-percentile. We summarize 
the forecast for observing a signal in neutrino telescopes or the scattering cross sections at 
future recoil experiments in Table [I] 

We have also presented model-independent parameterizations of the SI and SD cross 
sections for DM of arbitrary spins. 
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TABLE I: Survey of models and DM scattering predictions organized by the intrinsic spin of the 
DM candidate. The forecast for IceCube detection is also given, wherein we describe in what 
portion of the parameter space a 5a discovery can be made. Instances where the likelihood of 
detection via IceCube are low are marked by x. The 2a ranges for a^j ~ p and &sd~ P are gi ven 
in pb units. 
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4 x 10" 
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-13 _ 


4 x 10" 9 
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APPENDIX A: NUCLEAR MATRIX AND FORM FACTOR 

The spin dependent (SD) cross section involves the matrix element of spatial components 
of the quark axial-current 97759. The matrix elements of the proton (p) and the neutron 
(n) are 

(p\QTYsq\p) = 2Ajf ) (p\s p \p) , (n\qj-f 5 q\n) = 2Aj n) (n\s n \n) , (Al) 



where A q me asure s the fraction of the spin (s q = ^97759) carried by the quark q in the 



nucleon as in 



140| |. Isospin symmetry relates and A q n \ i.e. = A^ etc. Fro m the 



Wigner-Eckardt theorem, the matrix element for a nucleus (N) has the structure [1411 ] 

(N\q ll5 q\N) = 2\W(N\J N \N) . (A2) 
For an isolated single proton or neutron, the coefficients A a are simply 

AJM = AW , = AW . (A3) 



The numerical values 1421 ] for a proton are 

Ajf } = 0.78 ± 0.04 , Ajf) = -0.48 ± 0.04 , A( p) = -0.15 ± 0.04 . (A4) 
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The Wigner-Eckardt coefficient X[ N ^ relates the quark spin matrix element to angular mo- 
mentum of the nucleons. For the case of a nucleus with a single unpaired valence-nucleon 
(p,n), the shell model gives the Lande formula, 

A, = \{l + W + 1) - £(£ + 1)]/[J( J + 1)]}A^) . (A5) 

For general cases the estimates of the coefficients are parametrized as 

Af^(l/J)[Af)(5 p )+AW(5 n )], (A6) 



as in the review of Ref. 



129| | . The quantities (S p ) and (S n ) represent the respective expec- 



tation values for the pr oton and neutron group spin contents in the nucleus. They are very 



model dependent. Ref. [1431 ] evaluate d (S p ) = 0.03, (S n ) = 0.378 for 73 Ge. There is also a 
multiplicative form factor correction 1431 ] F(Q 2 ) to the scattering amplitude, that can be 



straightforwardly incorporated. The momentum transfer squared is 

Q 2 = \Pn ~p'n\ 2 = m l v l sin2 ^ > ( A7 ) 

with m r the reduced mass, v r the relative velocity of the DM particle and the target nucleus 
and 9 the scattering angle. It is expected that the form factor effects are only relevant for a 
sizable nucleus. Usually, a strong isospin decomposition is assumed in the parameterization. 
The isoscalar and isovector parts are 

a = (A u + A d )(X u + A d ) + 2A S X S , a ± = ( A u - A d ) (A u - X d ) . (A8) 

S(Q 2 ) = a 2 S 00 (Q 2 ) + a 2 S n (Q 2 ) + aoa^Q 2 ) (A9) 

Those S functions are well documented in nuclear model analyses. The form factor squared 
is given by 

|F(g 2 )| 2 = ^(g 2 )/,s(o) (aio) 

APPENDIX B: SPIN-1 AMPLITUDE OF B(k, fx) + f(p) -> B(£,u)f(q) 

We use the mUED model to illustrate the structure of the scattering amplitude. The 
LHT model follows a similar pattern. 
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The DM vector field Bi of odd parity couples to the known fermion / and its KK parity- 
odd partner F in the vertex BfF. The Feynman amplitude for the process B(k, fj,) + f(p) — > 
B(£,v)f(q) is 

- Cu(q)Y j-^- ^~z^2 ^ u (p) + ( crossing (i++u ,h++ -£) (Bl) 

where \i and v are the polarization indices. The heavy fermion F of mass M F is exchanged 
in the s channel (the first term), and in the u channel (the following parenthesis). We take 
the heavy M F 3> p, q, £, k limit for which 

(C/M F )u(q) A)Y + l»{jp- l)l u + 2M F gT) «(p) (B2) 

The 2M F term flips chirality. If the KK mode interaction is chiral, just like the SM, the term 
2M F vanishes when chirality projections are added in the amplitude. We use the Chisholm 
identity, 

VW = <TV ~ f v l a + 9 av Y + i^lfH* , (B3) 
and add the s and u channels contributions, 

(C/M 2 F )u(q) ((2p + k - ty^ + {2p + k- tYl v - (2 i>+ Jk- figT + i^ v \h + £) aW ) u{p) 

(B4) 

Then using the 4- momentum conservation, the Dirac equation of the massless quark limit, 
and the physical polarization properties, we obtain 

(C/M 2 F )u(q) ((p + g) V + (p + qy-f + ie^(k + £) aW ) u{p) . (B5) 

The relevant operators in the chirality basis are 

B»B V [q i(d V + d V)l(l ± Ts)?] and ie^{B v id a B^) q 1/3 (l ± l5 )q . (B6) 

The SD interaction is given by 

O sd = ie^(BJd a B,) qrtfnzq . (B7) 



[1] H. Baer and X. Tata, (2008), arxiv:0805.1905. 

[2] D. Hooper and E. A. Baltz, (2008), arxiv:0802.0702. 

34 



[3] S. Arrenberg, L. Baudis, K. Kong, K. T. Matchev, and J. Yoo, (2008), 0805.4210. 

[4] Eighth UCLA Symposium: Sources and Detection of Dark Matter and Dark Energy in the 

Universe , http://www.physics.ucla.edu/hep/dm08/. 
[5] 2nd International Workshop on the Interconnection Between Particle Physics and Cosmology 

, http://ppc08.phys.unm.edu/. 
[6] E. Komatsu et al, (2008), arxiv:0803.0547. 

[7] C. Boehm and P. Fayet, Nucl. Phys. B683, 219 (2004), hep-ph/0305261. 
[8] P. Fayet, Phys. Rev. D70, 023514 (2004), hep-ph/0403226. 
[9] M. Adibzadeh and P. Q. Hung, (2008), 0801.4895. 
[10] I. Antoniadis, Phys. Lett. B246, 377 (1990). 

[11] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999), hep-ph/9905221. 
[12] T. Appelquist, H.-C. Cheng, and B. A. Dobrescu, Phys. Rev. D64, 035002 (2001), hep- 
ph/0012100. 

[13] N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, Phys. Rev. D59, 086004 (1999), hep- 
ph/9807344. 

[14] P. Nath, Y. Yamada, and M. Yamaguchi, Phys. Lett. B466, 100 (1999), hep-ph/9905415. 
[15] D. A. Dicus, C. D. McMullen, and S. Nandi, Phys. Rev. D65, 076007 (2002), hep- 
ph/0012259. 

[16] G. Servant and T. M. P. Tait, Nucl. Phys. B650, 391 (2003), hep-ph/0206071. 
[17] C. Macesanu, Int. J. Mod. Phys. A21, 2259 (2006), hep-ph/0510418. 
[18] K. Kong and K. T. Matchev, JHEP 01, 038 (2006), hep-ph/0509119. 
[19] D. Hooper and S. Profumo, Phys. Rept. 453, 29 (2007), hep-ph/0701197. 
[20] I. Low, JHEP 10, 067 (2004), hep-ph/0409025. 

[21] J. Hubisz and P. Meade, Phys. Rev. D71, 035016 (2005), hep-ph/0411264. 

[22] V. Silveira and A. Zee, Phys. Lett. B161, 136 (1985). 

[23] J. McDonald, Phys. Rev. D50, 3637 (1994), hep-ph/0702143. 

[24] C. P. Burgess, M. Pospelov, and T. ter Veldhuis, Nucl. Phys. B619, 709 (2001), hep- 
ph/0011335. 

[25] V. Barger, P. Langacker, M. McCaskey, M. J. Ramsey-Musolf, and G. Shaughnessy, Phys. 

Rev. D77, 035005 (2008), arxiv:0706.4311. 
[26] H. Davoudiasl, R. Kitano, T. Li, and H. Murayama, Phys. Lett. B609, 117 (2005), hep- 



35 



ph/0405097. 

[27] B. A. Dobrescu, D. Hooper, K. Kong, and R. Mahbubani, JCAP 0710, 012 (2007), 
arxiv:0706.3409. 

[28] E. M. Dolle and S. Su, Phys. Rev. D77, 075013 (2008), 0712.1234. 

[29] S. Chang, L. J. Hall, and N. Weiner, Phys. Rev. D75, 035009 (2007), hep-ph/0604076. 

[30] D. Tucker-Smith and N. Weiner, Phys. Rev. D72, 063509 (2005), hep-ph/0402065. 

[31] K. Ishiwata, S. Matsumoto, and T. Moroi, (2007), 0710.3624. 

[32] T. Han and R. Hempfling, Phys. Lett. B415, 161 (1997), hep-ph/9708264. 

[33] J. R. Ellis, K. A. Olive, Y. Santoso, and V. C. Spanos, Phys. Rev. D71, 095007 (2005), 
hep-ph/0502001. 

[34] J. Ellis, K. A. Olive, and C. Savage, Phys. Rev. D77, 065026 (2008), 0801.3656. 

[35] CDMS, Z. Ahmed et al, (2008), arXiv:0802.3530. 

[36] J. Angle et al, Phys. Rev. Lett. 100, 021303 (2008), arxiv:0706.0039. 

[37] E. Aprile, 2008, presentation at UCLA DM 2008. 

[38] H. Wang, 2008, presentation at UCLA DM 2008. 

[39] J. Nikkei, 2008, presentation at UCLA DM 2008. 

[40] WARP, R. Brunetti et al, (2004), http://warp.lngs.infn.it. 

[41] ArDM, M. Laffranchi and A. Rubbia, (2007), hep-ph/0702080. 

[42] ZEPLIN-II, G. J. Alner et al, Phys. Lett. B653, 161 (2007), arXiv:0708.1883. 

[43] UK Dark Matter, G. J. Alner et al, Phys. Lett. B616, 17 (2005), hep-ex/0504031. 

[44] COUPP, W. J. Bolte et al, J. Phys. Conf. Ser. 39, 126 (2006). 

[45] KIMS, H. S. Lee. et al, Phys. Rev. Lett. 99, 091301 (2007), 0704.0423. 

[46] Super-Kamiokande, A. Habig, (2001), hep-ex/0106024. 

[47] Super-Kamiokande, S. Desai et al, Phys. Rev. D70, 083523 (2004), hep-ex/0404025. 

[48] A. Manalaysay, 2008, presentation at UCLA DM 2008. 

[49] PICASSO, M. Barnabe-Heider et al, Phys. Lett. B624, 186 (2005), hep-ex/0502028. 

[50] COUPP, E. Behnke et al, Science 319, 933 (2008), 0804.2886. 

[51] V. D. Barger, F. Halzen, D. Hooper, and C. Kao, Phys. Rev. D65, 075022 (2002), hep- 
ph/0105182. 

[52] M. Cirelli et al, Nucl. Phys. B727, 99 (2005), hep-ph/0506298. 

[53] F. Halzen and D. Hooper, Phys. Rev. D73, 123507 (2006), hep-ph/0510048. 

36 



[54] M. C. Gonzalez-Garcia, F. Halzen, and M. Maltoni, Phys. Rev. D71, 093010 (2005), hep- 
ph/0502223. 

[55] V. D. Barger, W.-Y. Keung, and G. Shaughnessy, (2007), arxiv:0709.3301. 

[56] V. Barger, W.-Y. Keung, G. Shaughnessy, and A. Tregre, Phys. Rev. D76, 095008 (2007), 

arXiv: 0708. 1325. 
[57] A. Gould, Astrophys. J. 388 (1992). 
[58] DAMA/LIBRA, R. Bernabei et al, (2008), 0804.2741. 
[59] B. Altunkaynak, M. Holmes, and B. D. Nelson, (2008), 0804.2899. 

[60] L. Roszkowski, R. Ruiz de Austri, and R. Trotta, JHEP 07, 075 (2007), arxiv:0705.2012. 
[61] E. A. Baltz, M. Battaglia, M. E. Peskin, and T. Wizansky, Phys. Rev. D74, 103521 (2006), 
hep-ph/0602187. 

[62] A. Raferty and S. Lewis, How many iterations in the gibbs sampler?, in Bayesian Statistics, 
edited by J. Bernardo, J. Berger, A. Dawid, and A. Smith Vol. 4, pp. 763-773, Oxford 
University Press, 1992. 

[63] A. Raferty and S. Lewis, Statistical Science 7, 493 (1992). 

[64] A. Raferty and S. Lewis, The number of iterations, convergence diagnostics and generic 

metropolis algorithms, in Practical Markov Chain Monte Carlo, edited by D. S. W.R. Gilks 

and S. Richardson, Chapman and Hall, 1995. 
[65] J. Dunkley, M. Bucher, P. G. Ferreira, K. Moodley, and C. Skordis, Mon. Not. Roy. Astron. 

Soc. 356, 925 (2005), astro-ph/0405462. 
[66] R. R. de Austri, R. Trotta, and L. Roszkowski, JHEP 05, 002 (2006), hep-ph/0602028. 
[67] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, (2008), arxiv:0803.2360. 
[68] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput. Phys. Commun. 149, 103 

(2002), hep-ph/0112278. 

[69] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput. Phys. Commun. 174, 577 

(2006) , hep-ph/0405253. 

[70] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput. Phys. Commun. 176, 367 

(2007) , hep-ph/0607059. 

[71] H. Baer and X. Tata, Weak scale super symmetry: From superfields to scattering events , 

Cambridge, UK: Univ. Pr. (2006) 537 p. 
[72] P. Binetruy, Super symmetry: Theory, experiment and cosmology , Oxford, UK: Oxford Univ. 

37 



Pr. (2006) 520 p. 

[73] M. Drees, R. Godbole, and P. Roy, Theory and phenomenology of sparticles: An account 
of four- dimensional N=l supersymmetry in high energy physics , Hackensack, USA: World 
Scientific (2004) 555 p. 

[74] H. Baer, C.-h. Chen, F. Paige, and X. Tata, Phys. Rev. D52, 2746 (1995), hep-ph/9503271. 
[75] K. L. Chan, U. Chattopadhyay, and P. Nath, Phys. Rev. D58, 096004 (1998), hep- 
ph/9710473. 

[76] J. L. Feng, K. T. Matchev, and T. Moroi, Phys. Rev. D61, 075005 (2000), hep-ph/9909334. 
[77] H. Baer, V. Barger, G. Shaughnessy, H. Summy, and L.-t. Wang, Phys. Rev. D75, 095010 
(2007), hep-ph/0703289. 

[78] N. Arkani-Hamed, A. Delgado, and G. F. Giudice, Nucl. Phys. B741, 108 (2006), hep- 
ph/0601041. 

[79] H. Baer, A. Mustafayev, E.-K. Park, and X. Tata, JCAP 0701, 017 (2007), hep-ph/0611387. 
[80] M. Drees and M. M. Nojiri, Phys. Rev. D47, 376 (1993), hep-ph/9207234. 
[81] H. Baer and M. Brhlik, Phys. Rev. D53, 597 (1996), hep-ph/9508321. 
[82] H. Baer and M. Brhlik, Phys. Rev. D57, 567 (1998), hep-ph/9706509. 
[83] P. Nath and R. Arnowitt, Phys. Rev. Lett. 70, 3696 (1993), hep-ph/9302318. 
[84] V. D. Barger and C. Kao, Phys. Rev. D57, 3131 (1998), hep-ph/9704403. 
[85] J. R. Ellis, T. Falk, and K. A. Olive, Phys. Lett. B444, 367 (1998), hep-ph/9810360. 
[86] J. R. Ellis, T. Falk, K. A. Olive, and M. Srednicki, Astropart. Phys. 13, 181 (2000), hep- 
ph/9905481. 

[87] J. R. Ellis, K. A. Olive, and Y. Santoso, Astropart. Phys. 18, 395 (2003), hep-ph/0112113. 
[88] J. R. Ellis, K. A. Olive, Y. Santoso, and V. C. Spanos, Phys. Lett. B565, 176 (2003), 
hep-ph/0303043. 

[89] H. Baer and C. Balazs, JCAP 0305, 006 (2003), hep-ph/0303114. 

[90] A. Djouadi, M. Drees, and J.-L. Kneur, JHEP 03, 033 (2006), hep-ph/0602001. 

[91] V. Barger, D. Marfatia, and A. Mustafayev, (2008), arxiv:0804.3601. 

[92] CDF/DO, (2007), hep-ex/0703034. 

[93] G. Bertone, D. G. Cerdeno, J. I. Collar, and B. C. Odom, Phys. Rev. Lett. 99, 151301 

(2007), arxiv:0705.2502. 
[94] B. S. Acharya, K. Bobkov, G. L. Kane, J. Shao, and P. Kumar, (2008), 0801.0478. 



38 



[95] C. Panagiotakopoulos and K. Tamvakis, Phys. Lett. B446, 224 (1999), hep-ph/9809475. 
[96] A. Dedes, C. Hugonie, S. Moretti, and K. Tamvakis, Phys. Rev. D63, 055009 (2001), hep- 
ph/0009125. 

[97] S. Y. Choi, D. J. Miller, and P. M. Zerwas, Nucl. Phys. B711, 83 (2005), hep-ph/0407209. 
[98] J. F. Gunion, D. Hooper, and B. McElrath, Phys. Rev. D73, 015011 (2006), hep-ph/0509024. 
[99] V. Barger, P. Langacker, H.-S. Lee, and G. Shaughnessy, Phys. Rev. D73, 115010 (2006), 
hep-ph/0603247. 

[100] V. Barger, P. Langacker, and G. Shaughnessy, Phys. Lett. B644, 361 (2007), hep- 
ph/0609068. 

[101] V. Barger, P. Langacker, and G. Shaughnessy, Phys. Rev. D75, 055013 (2007), hep- 
ph/0611239. 

[102] E. Accomando et al, (2006), hep-ph/0608079. 

[103] F. Ferrer, L. M. Krauss, and S. Profumo, Phys. Rev. D74, 115007 (2006), hep-ph/0609257. 
[104] V. Barger et al, Phys. Rev. D75, 115002 (2007), hep-ph/0702036. 

[105] C. Balazs, M. S. Carena, A. Freitas, and C. E. M. Wagner, JHEP 06, 066 (2007), 
arxiv:0705.0431. 

[106] H.-S. Lee, K. T. Matchev, and T. T. Wang, Phys. Rev. D77, 015016 (2008), 0709.0763. 
[107] U. Ellwanger, J. F. Gunion, C. Hugonie, and S. Moretti, (2004), hep-ph/0401228. 
[108] S. Chang, P. J. Fox, and N. Weiner, JHEP 08, 068 (2006), hep-ph/0511250. 
[109] V. Barger, P. Langacker, M. McCaskey, M. J. Ramsey-Musolf, and G. Shaughnessy, in 
preparation. 

[110] J. McDonald, N. Sahu, and U. Sarkar, (2007), arxiv:071 1.4820. 

[Ill] S. Profumo, M. J. Ramsey-Musolf, and G. Shaughnessy, JHEP 08, 010 (2007), 
arxiv: 0705. 2425. 

[112] G. Belanger, A. Pukhov, and G. Servant, JCAP 0801, 009 (2008), 0706.0526. 
[113] D. Hooper and G. Servant, Astropart. Phys. 24, 231 (2005), hep-ph/0502247. 
[114] N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, Phys. Lett. B429, 263 (1998), hep- 
ph/9803315. 

[115] T. G. Rizzo and J. D. Wells, Phys. Rev. D61, 016007 (2000), hep-ph/9906234. 

[116] I. Gogoladze and C. Macesanu, Phys. Rev. D74, 093012 (2006), hep-ph/0605207. 

[117] N. Arkani-Hamed, A. G. Cohen, and H. Georgi, Phys. Lett. B513, 232 (2001), hep- 

39 



ph/0105239. 

[118] T. Han, H. E. Logan, B. McElrath, and L.-T. Wang, Phys. Rev. D67, 095004 (2003), 
hep-ph/0301040. 

119] M. Perelstein, Prog. Part. Nucl. Phys. 58, 247 (2007), hep-ph/0512128. 
120] M. Schmaltz and D. Tucker- Smith, Ann. Rev. Nucl. Part. Sci. 55, 229 (2005), hep- 
ph/0502182. 

121] H.-C. Cheng and I. Low, JHEP 08, 061 (2004), hep-ph/0405243. 
122] H.-C. Cheng and I. Low, JHEP 09, 051 (2003), hep-ph/0308199. 

123] A. Belyaev, C.-R. Chen, K. Tobe, and C. P. Yuan, Phys. Rev. D74, 115020 (2006), hep- 
ph/0609179. 

124] C. T. Hill and R. J. Hill, Phys. Rev. D76, 115014 (2007), arxiv:0705.0697. 

125] V. Barger, W.-Y. Keung, and Y. Gao, Phys. Lett. B655, 228 (2007), 0707.3648. 

126] J. Hubisz, P. Meade, A. Noble, and M. Perelstein, JHEP 01, 135 (2006), hep-ph/0506042. 

127] Y. Bai, (2008), 0801.1662. 

128] G. Jungman and M. Kamionkowski, Phys. Rev. D51, 328 (1995), hep-ph/9407351. 
129] G. Jungman, M. Kamionkowski, and K. Griest, Phys. Rept. 267, 195 (1996), hep- 
ph/9506380. 

130] IceCube, A. Achterberg et ai, Phys. Rev. D76, 027101 (2007), arxiv: 0705. 1781. 

131] J. R. Ellis, A. Ferstl, and K. A. Olive, Phys. Lett. B481, 304 (2000), hep-ph/0001005. 

132] K. Cheung and T.-C. Yuan, JHEP 03, 120 (2007), hep-ph/0701107. 

133] A. D. Dolgov and S. H. Hansen, Astropart. Phys. 16, 339 (2002), hep-ph/0009083. 

134] J. R. Ellis and R. A. Flores, Nucl. Phys. B307, 883 (1988). 

135] V. Barger, W.-Y. Keung, H. E. Logan, and G. Shaughnessy, Phys. Rev. D74, 075005 (2006), 

hep-ph/0608215. 
136] K. S. Babu and E. Ma, (2007), arXiv:0708.3790. 

137] N. Kitazawa, N. Okada, and S. Sasaki, Phys. Lett. B380, 324 (1996), hep-ph/9602204. 
138] N. Arkani-Hamed, A. G. Cohen, E. Katz, and A. E. Nelson, JHEP 07, 034 (2002), hep- 
ph/0206021. 

139] J. Hubisz, J. Lykken, M. Pierini, and M. Spiropulu, (2008), 0805.2398. 
140] J. R. Ellis, R. A. Flores, and S. Ritz, Phys. Lett. B198, 393 (1987). 
141] M. W. Goodman and E. Witten, Phys. Rev. D31, 3059 (1985). 

10 



[142] G. K. Mallot, (1999), hep-ex/9912040. 

[143] V. Dimitrov, J. Engel, and S. Pittel, Phys. Rev. D51, 291 (1995), hep-ph/9408246. 



41 



